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Abstract. In this proceeding, recent theoretical investigations by the authors on the multiferroic _RMn03 
perovskites are briefly reviewed at first. Using the double-exchange model, the realistic spiral spin order 
in undoped manganites such as TbMn03 and DyMnOs is well reproduced by incorporating a weak next- 
nearest neighbor superexchange (~ 10% of nearest neighbor superexchange) and moderate Jahn- Teller 
distortion. The phase transitions from the A-type antiferromagnet (as in LaMnOa), to the spiral phase 
(as in TbMnOs), and finally to the E-type antiferromagnet (as in HoMnOa), with decreasing size of the 
R ions, were also explained. Moreover, new results of phase diagram of the three-dimensional lattice are 
also included. The ferromagnetic tendency recently discovered in the LaMnOa and TbMnOs thin films is 
explained by considering the substrate stress. Finally, the relationship between our double-exchange model 
and a previously used J1-J2-J3 model is further discussed from the perspective of spin wave excitations. 

PACS. 75.80.+q Magnetomechanical and magnetoelectric effects, magnetostriction - 75.47.Lx Manganites 
- 75.30.Kz Magnetic phase boundaries (including magnetic transitions, metamagnetism, etc.) - 75.30.Ds 
Spin waves 



1 Introduction 

Recently, the multiferroic materials, in which the ferro- 
electric (FE) and magnetic orders coexist and are inti- 
mately coupled, have attracted much attention due to 
their technological relevance and fundamental science chal- 
lenges [UHJIHlll]- Among the single phase multiferroic ma- 
terials, the undoped manganites with small size R (rare- 
earth) cations (like Tb, Dy, Ho and so on) is one of the 
most fascinating families. Not only the perovskites or hexag- 
onal i?Mn0 3 , BiMnOa, and YMn0 3 [31|Bl|Tl|gp[Ttj] . but 
also the i?Mn 2 Os series [TTlfT2lfT3] show multiferroicity. 

In this work, we only consider the undoped perovskite- 
type i?Mn0 3 . With decreasing R size, the ground state of 
these i?Mn0 3 compounds changes from the A-type anti- 
ferromagnet (A-AFM) (like LaMn0 3 and NdMn0 3 ), to 
the spiral spin state (like TbMn0 3 and DyMn0 3 ), and 
then finally to the E-AFM (like HoMn0 3 ) [6lfL4]. which 
is referred below as the "A-S-E transition" . On one hand, 
the spiral spin order (SSO) can break the space inver- 
sion symmetry and thus induce the observed ferroelectric 
(FE) polarization [51[f51[71[5Pl[T5] . although its microscopic 
mechanism remains under debate [16lfL7lfT8l[19] . This SSO 
driven improper ferroelectricity is also observed in other 
materials [20] • On the other hand, the E-AFM spin or- 
der can also induce the FE polarization because its zigzag 



chains break the space inversion symmetry as well [10 1 I2T 1 

mm- 

In contrast to the E-AFM phase which can be eas- 
ily obtained using the two-orbital double-exchange (DE) 
model |24j . the origin of the SSO remains a puzzle. A di- 
rect but phenomenological route to generate a SSO phase 
is via the magnetic frustration between NN ferromagnetic 
(FM) and next-nearest-neighbor (NNN) antiferromagnetic 
(AFM) interactions [3], e.g. via a J1-J2-J3 model with 
classical spins, where J\ is the NN superexchange (SE) 
while J 2 ( J 3 ) is the NNN SE along the b (a) directions [25] . 
Another route to obtain the SSO phase is to incorporate 
the Dzyaloshinskii-Moriya (DM) interaction (cx S,; x Sj) 
into the DE framework [TTII21)] . However, these models are 
not sufficient to describe the several phases of i?Mn0 3 . In- 
stead, our recent work has proposed an alternative model 
to understand the SSO and A-S-E transition in i?Mn0 3 
[27], which will be the focus of this manuscript. 

The rest of paper is organized as follow: In Sec. 2, 
we introduce the DE model and calculation methods. In 
Sec. 3, our recent study of the SSO and phase diagram in 
the two-dimensional (2D) lattice is briefly reviewed |27j . 
In Sec. 4, the study of the phase diagram is extended to 
the three-dimensional (3D) lattice. In Sec. 5, discussions 
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Fig. 1. (Color online) (a) Sketch of the crystal structure (a — b 
plane) of _RMn03. Two types of coordinate axes (a— b and x—y) 
are shown, (b) Illustration of the two kinds of distortions dis- 
cussed in the text: GdFe03-type (oxygen moves perpendicular 
to the Mn-O-Mn bond) and Jahn- Teller type (oxygen moves 
along Mn-O-Mn bond), (c) Two Jahn- Teller distortion modes: 
Q2 and Q3. Reproduced from [27] . Copyright 2008, the Amer- 
ican Physical Society. 

regarding the spin- wave spectrum of SSO are presented. 
The main conclusions are summarized in Sec. 6. 



2 Model and Methods 

Theoretically, the physics of manganites can be qualita- 
tively understood within the framework of the DE model, 
including the SE between the t 2g spins, and the Jahn- 
Teller (JT) interaction [55]. To obtain the spiral spin or- 
der, a weak NNN SE also has to be considered. The whole 
Hamiltonian reads as: 

a/3 

H = - Yl t* 13 Q iA<x C 3P + J AF Si ' S J 

<ij> <ij> 

+ ^ </2 7 Si • Sfc + \^2 / [Q2,iTx.i + Qz,iT z ,i] 
[ik] i 

+5 E^t,* + (!) 

i 

where the first term is the standard DE of the e g electrons. 
The DE hopping amplitudes tp are orbital- and direction- 
dependent. In particular, f*' 1 = t\j x — it 2 ' 2 — 3t 2 ' 2 = |t , 

j-1,2 _ j-2.1 _ +1,2 _ +2,1 _ \/3o xl,l _ 0.1,2 _ o2,l _ n 
b y — L y — b x ~ b x — 4 l 0j L z ~ °z ~ °z ~ U J 

and t 2 ' 2 = to where the superscript 1 (2) denotes the e g 
orbital d x 2 -y 3 (rf 3z 2_ r 2), and to (~ 0.2 — 0.3 eV) is taken 
as the energy unit [28]. The infinite Hund coupling used 
here generates a Berry phase fiy — cos(#i/2) cos(#,-/2) + 
sin(0j/2) sin(dj/2) exp[— i(tpi — <pj)], where 9 and ip are the 
angles defining the i 2g spins S in spherical coordinates. 
The second term in the Hamiltonian is the usual AFM 
SE coupling between NN i 2g spins. A realistic Jaf value 
is about ~ 10%to according to previous studies [IS]. The 



third term is the NNN SE between 3c? spins, where 7 is 
the direction index. J 2 is anisotropic due to the lattice 
distortion, e.g. the NNN SE J 2 b coupling between Mn(l) 
and Mn(4) can be ~ 1.8 — 2.2 times the value of the J 2a 
coupling between Mn(2) and Mn(3) since b > a [57], as 
sketched in Fig. 1(a). The fourth term is the electronic- 
phonon coupling of JT distortion, where A is the spin- 
phonon coupling coefficient and r is the orbital pseudospin 
operator, given by t x = c\c b + c\c a and t z =c^c a - c\c b 
[2"8j . For all i?Mn03 materials at low temperatures, |Q 2 | 
and Q3 are uniform with |Q 2 | w — The sign of Q2 

is staggered, which gives rise to the well-known staggered 
d 3x 2_ 1 .2 and d 3y 2_ r 2 orbital ordering |29j . The last term is 
the elastic energy of JT phonons. 

The numerical method used in this manuscript, and in 
the cited related references by our group, is the variational 
method employed at zero temperature (zero-T). The total 
energies (per site) of several candidate phases are com- 
pared to determine which is the most likely ground state. 
The 2D candidate phases include: A- AFM, C-AFM, CE, 
C1/4E3/4, Ci /3 E 2/3 , Dimer, E-AFM, G-AFM, SSO state 
24 , 27 , 28 1 . In this list almost all the 2D typical spin or- 
der patterns discussed in manganites have been included. 
For the SSO, spirals with wave vectors q from to 1/4 
are taken into account [30]. For the 3D calculation, the 
FM phase and canting spin state are added to the can- 
didate list. In principle, the t 2g spins S in Eq.(l) can be 
Heisenberg-like. However, in real manganites, all spin pat- 
terns of the known magnetically ordered phases can be de- 
scribed using the X-Y model, namely considering the exis- 
tence of an easy magnetic plane, such as the b-c plane for 
the SSO in TbMn0 3 and DyMn0 3 . Therefore, in practice, 
all the candidate phases considered here have X-Y model 
like spin patterns. Once a spin pattern and a JT distortion 
are selected, the total energy (per site) can be calculated 
in the infinite size lattice limit (thus, there are no finite- 
size effects here). This energy includes: (1) the DE kinetic 
energy Ek (including also the JT contribution) obtained 
from the exact diagonalization of the first and fourth terms 
in Eq.(l); (2) the SE energy Ej directly calculated from 
the second and third terms in Eq.(l). The elastic energy of 
the JT lattice distortions will not be taken into consider- 
ation since the Q 2 and Q 3 degrees of freedom are fixed in 
our variational method (therefore the last term of Eq.(l) 
is just the same constant for all candidate phases). For 
more details of the Hamiltonian and numerical methods, 
readers should consult Ref. [28] . 



3 Two-dimensional results 

In our recent publication [27] . the 2D phase diagram was 
studied, explaining the origin of the realistic spiral order 
and A-S-E transition in i?Mn0 3 . The main results can be 
summarized as follows: 

(1) With only the DE plus NN SE interactions (the 
lst-2nd terms of Hamiltonian Eq. 1), there is no spiral 
phase existing between the A- AFM and E-AFM, agreeing 
with previous studies [24] , 
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(2) By considering the DE, NN SE and NNN SE inter- 
actions (the lst-3th terms of Hamiltonian Eq. 1), the SSO 
phase emerges between the A-AFM and E-AFM when the 
NNN SE J 2b is larger than 0.017i , as shown in Fig. 2(a). 
Therefore, the A-S-E transition can be qualitatively un- 
derstood as the enhancement of Jaf /to and J2A0 by the 
GdFeC>3 distortion. However, the required J2 for the SSO 
remains too large and the obtained q for SSO is lower than 
the real value in i?Mn03. 

(3) To reproduce the realistic SSO, the whole Hamil- 
tonian (Eq. 1) should be considered, including the contri- 
bution from JT distortions. With a modest JT distortion 
coupling (AIQ2I = 1-5), the SSO phase region is expanded 
in parameter space, as shown in Fig. 2(b-c) [27] ■ The re- 
alistic short-wavelength SSO in i?Mn03 can be obtained 
with a weak J 2 b (~ 10% Jaf ~ l%io)- In addition, the JT 
distortion contributes to the insulating nature of i?Mn03 , 
as show in Fig. 2(d), which is crucial for the FE polariza- 
tion. 

(4) The phase transition between the A-AFM and SSO 
phases is second-order because the wave vector q changes 
continuously from (A-AFM) to a finite value (SSO). In 
contrast, the S-E and A-E phase transitions are of first- 
order. These orders of the phase transitions are indepen- 
dent of the JT distortions. Interesting physical phenom- 
ena, such as the bicritical point and phase separation, may 
emerge in the vicinity of the S-E phase boundary. 

(5) The Monte Carlo (MC) simulation on a 12 x 12 
lattice confirms the stability of the spiral phase and A- 
S-E transition at low temperature. The temperature de- 
pendent FE polarization is also obtained using the phc- 
nomenological equation — e,-^ x (Si x Sj), which agrees 
with experimental observation qualitatively. Note here that 
a finite-size lattice is used in our MC simulations, and 
those finite lattices can only accommodate some particu- 
lar SSO, such as the state with q = 1/6 wavevector (cor- 
responding to Tbo.4iDyo.5gMn03 [9]). In particular, the 
12x12 lattice is a very good choice since it is compatible 
with the A-AFM and E-AFM states, as well as the q = 1 /6 
SSO simultaneously. At present, MC simulation results on 
larger 2D or 3D lattices is not available due to the rapid 
growth of CPU time with increasing lattice sizes. There- 
fore, in this Proceeding, we will focus most of our attention 
on the ground states using the zero-T variational method. 
Readers can consult our original publication (Ref. [27]) if 
they are interested in the finite-temperature results on the 
2D 12 x 12 lattice. 

Note that this weak NNN SE interaction (< 10% NN 
SE), while shown here to be crucial in the context of the 
manganite multiferroics, it does not alter the previous 
large body of investigations and conclusions reached via 
MC simulations for undoped and doped LaMn03, since in 
that case the extra NNN SE couplings can be neglected. 

4 Three-dimensional results 

In the above 2D study, the ferromagnetic (FM) phase was 
neglected in the candidate list. This was reasonable since 
none of the i?Mn03 bulk materials is FM. However, very 
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Fig. 2. (Color online) (a-c) Zero-T phase diagram of the 2D 
two-orbital DE model for 7?Mn0 3 . Notations: A: A-AFM; C: 
C-AFM; D: Dimer; E: E-AFM; G: G-AFM; S: Spiral spin state, 
(a) Results without the JT distortion. The possible path for the 
A-S-E phase transition is indicated by the arrow. The phase 
diagram is independent of as long as J2& < Jih- (b) With 
the JT phonon (A|Q 2 | = 1.5) and J 2a = 0. The q's for TbMn0 3 
and DyMnOs are also indicted by broken lines, (c) The same 
with (b) except J 2a ~ 0.5 J2&. In (a-c), all the slopes of the A-S 
boundaries are —2 since the SE energy Ej for the A-AFM and 
SSO phases is 2Jaf cos(5) + J 2 h cos(2<5) + J2a, where 5 is the NN 
spin angle, (d) Density of states for the A-AFM, SSO (q = 0.14 
as TbMn0 3 ), and E-AFM with the same JT distortion, (a-c) 
Reproduced from [27]; Copyright 2008, the American Physical 
Society. 

recently, the FM tendency in LaMnOs and TbMn03 thin 
films grown on SrTiOs substrate has been observed [311 
32,33]. To fully understand this exotic behavior, a 3D lat- 
tice is necessary. In addition to the candidate phases in 
the above 2D studies, two more phases: the FM and 3D 
canting spin phase are included. The 3D canting spin or- 
der is similar to the forementioned SSO, namely both of 
them are noncollinear spin orders and can be scaled by 
wave vectors q. However, the SSO is only noncollinear in 
the a — b plane, while spins are collinear along c axis. In 
contrast, in the 3D canting spin order, spins are noncolin- 
ear along all the directions, namely the NN spins' angles 
are isotropic. To simplify the model, all the 12 NNN SE 
couplings are firstly assumed to be equal. 

Using the same zero-T variational method, we first 
calculate the phase diagram on an infinite cubic lattice 
without the JT distortion, as shown in Fig. 3(a). Com- 
paring with Fig. 2(a), there are several differences. In 
the 3D lattice phase diagram, the FM phase occupies a 
large region at low Jaf and J 2- The A-AFM can not ex- 
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ist until J 2 > 0.018. This FM-rich phase diagram can 
shed light to the FM tendency in i?MnC>3 thin films. The 
LaMnC>3 lattice on a SrTi03 substrate is almost cubic 31, 
[34] . Even for the TbMnC>3 thin film, the differences be- 
tween a, b, and c/y/2 are also reduced by the stress [321133] . 
namely the thin-film lattice becomes closer to cubic than 
that of bulk. This FM enhancement in the cubic lattice 
can be understood based on three observations: (1) the 
isotropic exchange interactions due to the isotropic Mn- 
O-Mn bond-length and bond-angle do not prefer phases 
with anisotropic spin patterns, e.g. A-AFM; (2) the Mn- 
O-Mn bond-angle is more straight in the higher-symmetry 
lattice, which will enhance the ratio between the DE and 
SE terms, giving rise to an increased bandwidth and FM 
tendency; (3) the orthorhombic distortion in the bulk is 
advantageous for the d 3x 2_ r 2 /d 3y 2_ r 2 type orbital-ordering 
in the A-AFM and SSO phases, but this orbital order will 
be suppressed in the nearly cubic lattice on the substrate. 
In summary, the FM tendency in i?Mn03 is natural once 
the lattice is close to cubic. 



Similar to the 2D case, to understand the properties of 
i?MnC>3 in the bulk, the JT distortion of orthorhombic lat- 
tice has to be considered. Using the same AIQ2I = 1-5, the 
phase diagram is recalculated, as shown in Fig. 3(b). The 
FM region is largely suppressed, and it is being replaced 
by a robust A-AFM phase, which agrees with the proper- 
ties of real i?MnC>3 (R —La, Pr, Nd, ...) bulk materials. 
The g's for TbMnC>3 is also indicted. However, the spiral 
phase region is somehow narrow, and the realistic g's for 
DyMnOs are missing. To solve this puzzle, the anisotropy 
of NNN SE has to be considered [27] . Since the accurate 
ratios between J2's along different directions are unclear, 
here only two limits are calculated. In Fig. 3(a-b), the 
upper limit, namely isotropic J2's, has been considered. 
With the isotropic J2, six equal NNN SE bonds per site 
are taken into account. For the lower limit, only the J2b 
is nonzero (it is along the 6-direction because it has the 
strongest intensity due to the largest Mn-O-O-Mn angle), 
while other J2S are all set to be zero. In this limit case, 
only one NNN SE bond per site is considered. The new 
phase diagram is calculated and shown in Fig. 3(c). The 
spiral phase region is expanded in this revised phase di- 
agram, and now the wave vector q for DyMnOs can be 
found in Fig. 3(c). In real manganites, the J 2 should be 
within these two limits, namely six inequivalence NNN SE 
bonds per site should be considered. 

In short, the 3D calculation agrees with the 2D results 
qualitatively. Besides the SSO and the A-S-E transition, 
the FM tendency in i?MnC>3 thin films is also explained by 
considering the substrate stress. The several phases exist- 
ing in the phase diagram illustrate the possibility to modu- 
late the subtle phase competition in i?Mn03 using various 
methods, such as stress and strain, or external magnetic 
and electric fields. With these stimulations, phase sepa- 
ration may be possible to emerge, which would result in 
colossal responses to these external stimulations. 
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Fig. 3. (Color online) Zero-T phase diagram of the 3D two- 
orbital DE model for i?Mn0 3 . Notations: F: FM; CT: 3D cant- 
ing spin order; others conventions are the same as those in 
Fig. 2(a-c). (a) Results without the JT distortion and with 
the isotropic J2 coupling, (b) Results with the JT distortion 
(AIQ2I = 1.5) and with the isotropic J2 coupling. The q's 
for TbMnOa are indicated, (c) Same as (b) except that the 
anisotropic J2 is used. Only the Jib is nonzero while others are 
zero. 

5 More Discussion 

In our publication [27J , we argued that although the clas- 
sical spin model with NNN magnetic frustration provides 
a natural starting point to describe phenomenologically 
the spiral phase, this simple model is not sufficient to 
deeply understand the microscopic origin of the SSO in 
perovskite manganites. In the J1-J2-J3 model with classi- 
cal spins, a strong J2 (compared with Ji) is required to 
generate a realistic short-wavelength SSO, e.g. to reach 
the realistic wave vectors of TbMnOs and DyMnOs, the 
h/\J\\ should be about 0.78 ~ 1.35 [25]. However, the 
first-principles calculations show that the value of J2 is 
only about 0.56 |Ji| [18], which is not sufficient to induce 
the experimentally observed short-length (L = 6) SSO. 
Also, the coupling J3 was found to be 0.375 |Ji| and in 
addition it is AFM, opposite to the weak FM coupling 
(O.OUi) used in the previous model study [25]. Further- 
more, the classical spin model is not suitable to explore 
in a single framework the several phases (A-AFM, spiral, 
and E-AFM) found in _RMn03, e.g. it can not distinguish 
the E-AFM order from the q = 1/4 spiral order since they 
have the same energy. 

Very recently, the electromagnon excitations in HMnOs 
were studied |35j|36j. The dynamics of the spiral order is 
a very important issue, which can be helpful to discover 
the real mechanism of multiferroicity. The J1-J2-J3 model 
with classical spins was used to calculate the spin wave 
properties [36j . Thus, it is necessary to clarify the proper- 
ties of the classical spin model further to decide whether 
such an approach is sufficient to study the spin waves in 
i?Mn0 3 . 
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Since the e g density is uniform for i?MnC>3, in the fol- 
lowing we can use a classical approximation for the DE 
interaction [37] : 

H = [Jde\/1 + Si ■ Sj + JafS.; • Sj] + J2 7 Si ■ Sj, 

<ij> [ij] 

(2) 

where Jde is an effective coupling which can be derived 
from the kinetic energy. Note that this approximation gives 
the same energy as in the full study of the DE model only 
when the NN spins' couplings are uniform. Also, this ap- 
proximation is accurate only when the JT distortions are 
neglected. For instance, Jde is about — 0.423io for the 2D 
A-AFM or the spiral order in i?MnC>3 since the DE kinetic 
energy (per site) is about — 1.196to- But for the 3D FM 
and canted states of i?MnC>3, Jde is about — 0.363io since 
the DE kinetic energy (per site) is about — 1.54U - In the 
following, we will only focus our attention on the 2D SSO 
phase. 

Using the standard expansion technique and rotation 
of the quantum projection axis, the linear spin-wave the- 
ory can be applied to the SSO phase [38]. It is straightfor- 
ward to obtain an equivalent NN coupling J\ — — , JpE 

2a/ 1+cos(<5) 

Jaf, where 5 is the (ground state) angle between NN 
the spins. This coupling Ji is weaker than to and Jaf, 
and comparable with J 2 . Therefore, the classical J1-J2-J3 
model appears suitable for this simple description to some 
extent. 

Even when the JT distortion is included, although the 
exact formula is unknown the classical DE term can still 
merge with the NN SE term, resulting in a weak effective 
Ji coupling. Therefore, it seems possible to use the J1-J2- 
J3 model with classical spins to study the spin wave in the 
SSO i?Mn03. However, once there are phases with nonuni- 
form NN spins couplings (e.g. E-AFM and CE phases), the 
DE Hamiltonian can not be simplified into a mere classi- 
cal Heisenberg model, thus several phases in manganites 
can not be accessed by using a pure spin model. The phys- 
ical reason is that the e g electrons are itinerant (due to 
the DE hopping) even in the undoped i?Mn03, while the 
pure spin model is based entirely on localized spins. This 
itinerant DE process is essential to understand the novel 
physics in manganites, for instance the charge/orbital or- 
dering in the zigzag chains of the E-AFM or CE phases 
[2"gll3"5] , and also the multiferroicity in the E-AFM state 
of HoMnOs [2Tj . Thus, it is not correct to investigate the 
existence of several competing phases in i?Mn03, such as 
the phase diagram and phase transitions, using just a pure 
spin classical approximation. 

6 Conclusion 

In conclusion, here we have provided a microscopic de- 
scription of the several competing spin orders in multi- 
ferroic i?Mn03 perovskites. The experimentally observed 
spiral order and FE transition can be obtained by incor- 
porating a weak NNN superexchange interaction and a 
Jahn- Teller distortion into the standard two-orbitals DE 



model for manganites. Several aspects of the experimen- 
tally known A-S-E phase transition with decreasing R size 
are well reproduced by including the GdFe03-type dis- 
tortion in our study. Furthermore, the FM tendency in 
i?Mn03 thin films is also explained by considering the 
substrate stress. The relationship between a previously 
studied classical spin model for multiferroics and our more 
fundamental model is also further discussed. 

This work was supported by the NSF (DMR-0706020) and 
the Division of Materials Science and Engineering, U.S. DOE, 
under contract with UT-Battelle, LLC. J.M.L. was supported 
by the National Key Projects for Basic Research of China 
(2006CB921802 and 2009CB929501) and the National Natu- 
ral Science Foundation of China (50832002). S. Y. was sup- 
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